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SUMMARY 

In a refrigerated seal the fluid to be sealed flows through a refrigerated 
housing or constriction. The fluid can be frozen to the housing during the 
transient phase. Under steady-state conditions the refrigerated seal proved to 
be a dynamic low-leakage seal. The concept was extended to pressure differ- 
ences of 6.9 MPa (1000 psl). 


INTRODUCTION 

One purpose of a seal Is to prevent fluid loss (l.e., to minimize 
leakage). However, In many Instances Improperly designed seals have caused 
failures and deratings of turbomachines. Thus a second and equally Important 
purpose of a seal Is the prevention of dynamic Instabilities (ref. 1). 

In many turbomachine applications, fluids with significantly different 
boiling points may exist In close proximity. Such fluids could provide the 
heat sink for a refrigerated seal as discussed In Hendricks (ref. 2*). Flow 
losses and dynamic parameters such as stiffness and damping can be actively 
controlled by either the formation of a solid or an Increase In the viscosity 
of the fluid within the narrow aperture. The rotating member must generate 
energy through viscous dissipation to balance the refrigeration load. 

In this study the self-sealing, high-shear flow passage approach of 
Hendricks (ref. 2) was extended to large pressure differences. 


*In a note, "Guidelines for Success In Freeze Plugging Pipe," C.L. 
Schelllng, Freeze Seal, Inc., New Castle, Delaware, describes a freezing method 
to repair valves in lines to 48 In. In diameter by freezing a plug In the line 
(Ref. Nuclear Safety, 1985). 



APPARATUS AND PROCEDURE 


An apparatus (fig. 1) was constructed to determine the effectiveness of a 
self-sealing system at large pressure differences. The fluid to be sealed was 
water and the refrigerant was nitrogen. Apertures of the small- and large- 
length-to-dlameter-ratlo (L/D) Borda types (ref. 3) were used as shaft seals 
and placed between bearing supports (fig. 2). With the shaft rotating, water 
was added to the reservoir, which leaked profusely. As the nitrogen cooled an 
annular region around each seal, the leakage stopped entirely as was also dem- 
onstrated In Hendricks (ref. 2) for an orifice and a Borda aperture. Less time 
was generally required to stabilize the rotor and to seal the system for the 
large-L/D Borda seals, as would be expected from their better heat transfer and 
larger contact area. For eccentric seals or those rubbed, misaligned, or badly 
damaged, the solid plug (Ice layer) built up quickly, and with sufficient 
refrigeration the leak stopped completely. 

Once the plug Is formed and maintained, the seal should not "wear" and 
should adjust to system perturbations unless the pressure difference becomes 
sufficient to fracture the plug. To Illustrate this, the water reservoir was 
filled to capacity, taking care to eliminate all trapped air, and valved off. 
Using gaseous nitrogen as the pressurant, the water reservoir pressure was 
Increased In 0.3-MPa ( 50— ps 1 ) Increments and valved off. If no pressure drop 
or leaks were observed, the pressure was Incremented. At 6.9 MPa (1000 ps 1 g) 
one section of the plug began to leak. 


ANALYSIS 
Simplified Model 

The analysis for a simplified model Is given In the appendix. It Is shown 
that In a transient thermal analysis with freezing the calculated clearance 
became 0.01 cm (0.004 In) and some leakage could be expected. The shear stress 
calculated from freeze-plug failure data for the small-L/D aperture, 1.03 MPa 
(150 psl), was In reasonable agreement with the ultimate stress as estimated 
from the behavior of concrete. 


Finite Element Model 

Although It was the small-L/D Borda aperture plug that failed, both long 
and short apertures were modeled (fig. 1). The finite element model (FEM) 
consisted of 882 and 182 elements for the large- and small-L/D apertures, 
respectively. In both cases, there were 14 elements In the radial direction, 

8 In the Ice, and 6 In the steel, leaving 63 and 13 axial elements, respec- 
tively (fig. 3). Between 0.715 and 0.750 cm (0.28 and 0.295 In) the steel 
tube was supported (l.e., fixed); the small-L/D tube started at 0.715 cm. The 
curvature of the cooling coll was not considered. 

The axial pressure profile through the seal was assumed to be linear. 

Thus loading was that of a thin Ice shell with a metal tube backing. 
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The strength of Ice was at best elusive (see appendix); however, from pre- 
vious experience, the material behavior of Ice was assumed to be similar to 
that of concrete. The modified Pauw relation as adopted by ACI** Is 

E c = 33 W Wf c (1) 

where E c Is the secant modulus of elasticity ( psl ) , W the unit weight 
(lb/ft 3 ), and f c the ultimate strength (psl). The allowable shear stress 
on a section for an ultimate strength design Is 

2<p f c < V cu < 3.5 <p f c = 3.5 <pE c /(33W 3 ^ 3 ) (2) 

For Ice, E c = lxlO 6 psl, <p = 0.85, W = 56 lb/ft 3 , and the shear stress range 
becomes 


122 psl < V cu < 214 psl 
0.84 MPa < v cu < 1 .48 MPa 


(3) 


RESULTS 

Numerical 

In both cases (large- and small-L/D plugs) the shear stress diminished 
rapidly In the axial direction from the high-pressure side to the low-pressure 
side of the plug seal with concentrations near the point of fixity (fig. 4). 

The shear stresses near the high-pressure side were 2.76 and 4.28 MPa (400 and 
620 psl) for the large- and small-L/D plug seals, respectively. The average 
shear stresses were approximately 0.69 and 1.38 MPa (100 and 200 psl), respec- 
tively. Anticipating the behavior of Ice to be similar to that of concrete, 
the previously estimated shear stresses were 0.84 to 1.48 MPa (122 to 214 psl), 
Indicating the small-L/D plug seal to be overstressed; both plug seals were 
overstressed near the high-pressure side. 


Experimental 

At 6.9 MPa (1000 pslg) the plug fractured and water leaked out. The 
fracture was In reasonable agreement with the simplified and finite element 
analyses. Instrumentation of the passage would have provided more details for 
modeling. 


Applications 

The applications and limitations of a refrigerated seal are discussed In 
Hendricks (ref. 2). Among the applications are precise seal apertures; quick 
fix methods; solid formation In hypersonic flows; redundant or emergency seals; 
enhanced stability; and secondary fluid Injection. The seals cannot be used 
as static seals unless auxiliary heat Is available and sufficient housing 


** American Concrete Institute Publication ACI 318-83. 
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(shaft) strength Is provided to restrain the volumetric expansion of Ice forma- 
tion. It may also be possible to adaptively control stiffness and damping 
coefficients by controlling refrigerant flow to minimize system vibrations 
(Including the first critical). 


SUMMARY OF RESULTS 

A proposed sealing concept where the fluid to be sealed flows through a 
refrigerated housing or constriction has been extended to pressure differences 
of 6.9 MPa (1000 psl). As the fluid cools. Its viscosity Increases and It 
becomes more resist to leakage. In some cases this may be sufficient to estab- 
lish the desired leakage rates, but In other applications the temperature must 
be lowered to below the freezing point of the fluid In the aperture In order 
to build and maintain a sealing plug. 

Under steady operating conditions the viscous dissipation provided by the 
rotating member maintains a thin liquid film between the shaft and the sealing 
plug. For the water-nitrogen system axial pressure differences to 6.9 MPa 
(1000 psl) were maintained. The concept provided a dynamic low-leakage Borda 
seal . 


A finite element method was used to determine the strains Induced In the 
sealing plug. The predictions were In reasonable agreement with the experi- 
mental data. 

Discussions with Cole (ref. 4) clearly Indicate the need for a thermo- 
mechanical creep analysis. 


SYMBOLS 

A area 

b aperture static clearance 

c aperture dynamic clearance 

C p specific heat 

D hydraulic diameter 

E c secant modulus of elasticity 
f c ultimate strength 

G mass flux 

H latent heat 

h heat transfer coefficient 

k thermal conductivity 
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L ratio of heat transported by refrigerant to that conducted through solid 
plug, Eq. (8) (also length) 

N ratio of heat Input from seal to that absorbed by refrigerant, Eq. (7) 

Nu Nusselt number, hD/k 

Pr Prandtl number, Cpn/k 

R shaft radius 

Re Reynolds number, GD/n 

S viscous dissipation 

T temperature 

t thickness 

V cu shear stress, ultimate design 

W weight per unit volume 

p density 

© dimensionless time In terms of convection, conduction, and latent 
energies at Interface, Eq . ( 6) 

e time 

n viscosity 

t stress 

Q angular velocity 

Subscripts : 
a coolant side 

amb ambient 

f fusion or solid 

1 seal fluid side 

m melting 

max maximum 

n normal 

s surface 
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xc cross section 

rz surface normal to r, direction parallel to z 

0 reference or stagnation 


APPENDIX - SIMPLIFIED ANALYSIS 
Thermal 

Consider a shaft seal where either the working fluid can solidify or a 
solid plug can be formed through secondary fluid Injection (fig. 5). Further, 
consider a seal configuration such that the clearance-to-radlus ratio 

c/R = (b - t)/R Is much less than 1, so that the problem can then be addressed 

as one-dimensional; refrigeration Is provided through a heat exchanger 
(fig. 6). In terms of this model the viscous dissipation becomes 

S = (4) 

For a fixed rotational speed G, S Is directly related to the viscosity 
n, which generally Increases as temperature decreases*, and Is Inversely pro- 
portional to the clearance or the amount of solid plug formed In the seal 

aperture. As an example, for a 25-mm (1 In) shaft turning at 18 000 rpm In 

water with cold gaseous nitrogen refrigerant: 

1. Cooling the aperture water from 38 to 10 °C (100 to 50 °F) nearly 
doubled the viscosity and Increased the axial pressure drop across the aperture 

2. Under equilibrium conditions for b - t of 0.01 mm, S became 6 W/cm 2 
(19 000 Btu/hr ft 2 ). Thus for some applications It may only be necessary to 
control the temperature of the fluid within the aperture to minimize leakage. 
Once a solid plug Is formed, the energy required to maintain the plug can be 
small. 


A model of the energy transport required to reduce the seal aperture fluid 
temperature and build a solid plug to provide a clearance of b - t Is Illus- 
trated In figure 7. An analysis for this geometry follows that presented by 
Rohsenow (ref. 6). The governing parametric relation Is 


© = In 


N - 1 

N(L + 1) - 1 



(5) 


^Although Infinite at the critical point, viscosity has a minimum In the 
thermodynamic critical region. At low temperatures (near the triple point) 
viscosity has a maximum for quantum fluids. The temperature effect Is usually 
strong for liquids, exp (a/T), and weak for gases, T 2/3 . Further, for a syn- 
thetic oil such as Mobil RL-714 Stock 509, a temperature change from 60 to 
20 °C (140 to 68 °F) Increases the viscosity by a factor of 4 (Braun (ref. 5)). 
One now must consider the stability of operating In the critical region with 
bearings, etc, as n has a minimum. 
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where & Is the dimensionless time to freeze the thickness b - t In terms of 
conduction to, convection from, and latent heat at the solid-fluid Interface. 


e = 


< T o- 

P Hk 


V 


e 


( 6 ) 


where N Is the ratio of the heat Input from the seal to that absorbed by the 
refrigeration, 


N = 


h i ( T o - ¥ 

h a < T f - V 


(7) 


and L Is the ratio of the heat transported by the refrigerant to that con- 
ducted through the solid 


Here p and k are properties of the solid plug and H Is the latent heat of 
fusion. The heat transfer coefficients h a and h^ were determined to 
sufficient accuracy by using the McAdams equation 

Nu = 0.023 Re°’ 8 Pr°- 4 (9) 

even though more sophisticated forms are available. 

For gaseous nitrogen refrigerant (T = -193 °C (-316 °F), P = 0.138 MPa 
(20 psla)) flowing at 3 m/s (10 ft/s) through a 4.6-mm (0.18-1n) diameter tube, 
the Reynolds number Re was 12 400 and h a was 70 W/m 2 K (12 Btu/hr ft °F). 

On the water side, for a leakage velocity of 1.8 m/s (6 ft/s) with a hydraulic 
diameter (2c) of 2.5 mm (0.1 In), the Reynolds number became 3500 and h^ was 
9000 W/m 2 K (1600 Btu/hr ft 2 °F). For large To - T f freezing will not occur. 
However, at Tq •* Tf, (a reasonable assumption for the Borda Inlet) an explicit 
solution Is found 


As an example, consider the apparatus of figure 5 (Hendricks (ref. 2)). 
Using equations (6), (8), and (10) and the properties of Ice, with b = 1.2 mm 
(0.05 In) (l.e., t = 0), an Ice plug should form In at least 1/2 min. The 
actual "pinch off" clearance b - t was self-limiting; for as soon as the 
clearance became very small, S became large and melted the Ice. An additional 
feature was that the maximum temperature occurred within the aperture channel 
and not at the rotating shaft, thereby preventing complete closure so long as 
the shaft was rotating. 


Mechanical 

The Ice plug can be represented as a cylinder of thickness t (clear- 
ance), length L (length of Borda), and radius R 0 (shaft radius), fully 
bonded to the small- or large-L/D Borda Inlet (fig. 8). Although the adhesive/ 
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cohesive strength can vary with temperature and surface characteristics, we 
assumed It to be uniform. One end of the cylinder was loaded at Pq and the 
other was at P am b* The axial pressure distribution In the clearance gap was 
parabolic, but previous work with fluid nitrogen and hydrogen has shown that a 
linear approximation Is satisfactory (Hendricks and Stetz (ref. 7)). As such 
the average uniform radial load Is (Pq + p amb)/ 2 * 

Under these conditions a concrete-footer-slab model was used to simulate 
the behavior of the Ice plug. Balancing the forces on the model gave 


tA s = ( Pq - p amb) A xc 


(ID 


and the shear at failure became 

t = = 1000 x = 1.03 MPa (150 psl) (12) 

The allowable shear based on ultimate strength, which In turn was based on 
modulus of elasticity, varied between 0.8 and 1.5 MPa (115 and 215 psl), which 
was In reasonable agreement with equation (12). 

Note that the pressure difference restrained by the Ice depended on seal 
length and clearance and on the compressive and adhesive/cohesive strength of 
the Ice. Although we assumed the Ice to be fully bonded to the steel, the 
strength depended on the rate of freezing and the temperature-pressure history 
as well as on substrate characteristics. 

Discussions with Cole (ref. 4) and his work demonstrate that grain size 
d, creep rate e, temperature T, applied load a, and confining pressure all 
affect the strength and failure of Ice (In our case type 1^). For example, at 
-5 °C (22 °F) the strength peaks near 10 MPa (1450 psl) for e 10~3 and 

decreases to 1 MPa for e 10~5 . At -196 °C (-320 °F), the strength range 

Is 21 to 48 MPa (3 to 7 ks 1 ) for P near ambient to over 160 MPa (23 ksl) for 
50 MPa < P < 150 MPa ( 7.2 ksl < P < 22 ksl), (Durham et al. (ref. 8)). When 
these values are applied to the seal, the solid near the water-ice Interface 
would be at high strain and low strength, and the solid near the nitrogen cool- 
ing coll would have a much higher strength. Increasing the confining pressure 
and surface area should significantly Increase the capacity of the Borda seal 
over that of the orifice seal. These discussions clearly Indicate the need for 
a thermomechanical analysis Including the effects of the parameters e, d, T, a, 
and P. 
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FLUID FILL LINE 



Fig. 1 Schematic of test apparatus for small- and large-Borda type of refrigerated seals. (Dimensions 
are in millimeters (inches). ) 







Fig. 3 Finite element model of freeze plug seal. (Dimensions are in centimeters (inches). ) 
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Fig. 4 Shear stress at ice-steel interface. 



Fig. 5 Schematic of experimental test apparatus for Borda and orifice types of refrigerated seals. (Dimensions are 
in millimeters (inches).) 
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STEADY-STATE MODEL OF SEAL APERTURE 

Fig. 6 Schematic of velocity and temperature profiles for refr igerated-seal 
configuration. 







SEAL FLUID 




r \\ \\ 

COOLANT 

Ta 

Ficj. 7 Temperature profile for refrigerated-seal model. 




Fig. 8 Slab model simulation of ice plug seal. 
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